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INTRODUCTION 


r 


Let  X and  Y be  two  random  variables  with  cumulative  distri- 
bution functions  F(x)  and  G(y)  respectively.  Suppose  Y is  the 
strength  of  a component  subject  to  a stress  X.  Then  the  component 
fails  if  at  any  moment  the  applied  stress  (or  load)  is  greater  than 
its  strength  or  resistance.  The  stress  is  a function  of  the  environ- 
ment to  which  the  component  is  subjected,  and  its  value  at  any  point 
of  time  is  considered  a random  variable.  The  strength  of  a component 
is  measured  by  the  stress  required  to  failure.  Strength  depends  on 
material  properties,  manufacturing  procedures  and  so  on.  If  the  com- 
ponents under  question  are'  mass  produced  and  their  selection  in  a 
given  system  is  assumed  to  be  made  at  random,  then  the  strength  should 
also  be  considered  a random  variable.  The  reliability  of  a component 
during  a given  period  [0,T]  is  taken  to  be  the  probability  that  its 
strength  exceeds  the  stress  during  the  entire  interval,  that  is,  the 
reliability  function  R is  given  by 

R = P(X  < Y) 

From  practical  considerations  it  is  desirable  to  draw  inference  about 
the  reliability  function.  The  problem  of  estimating  R has  been  con- 
sidered by  many  using  nonparamctric , Bayesian  and  parametric  approach. 
We  shall  present  a survey  of  available  results  and  consider  some  new 
results . 

The  above  model  was  first  considered  by  Birnbaum  (1956)  and  has 
since  found  an  increasing  number  of  applications  in  many  different 


ll 
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areas,  especially  in  the  structural  and  aircraft  industries. 

As  an  example,  consider  the  following  problem  discussed  by 
Lloyd  and  Lipow  (1962).  A solid  propellant  rocket  engine  is  succus 
fully  fired  provided  the  chamber  pressure  (X)  generated  by  igint 
stays  below  the  burst  pressure  (Y)  of  the  rocket  chamber.  If 
X <i  Y,  the  engine  blows  up  and  the  operation  is  a failure. 

Note  the  problem  of  inference  about  R = P(X  < Y)  is  similar  tc 
the  problem  of  estimation  of  P = P (X  > Y) , the  probability  of  fail- 
ure. So  one  can  either  talk  of  R,  or  of  P. 

2 . Nonparametric  approach 

Let  (X^,  X2,...,Xjj^)  and  (Yj^,  Y2,...,Y^)  be  two  independent 


samples 

of  measurements  on 

X and 

Y respectively 

. Let 

r 

- 

1 ACCrsSION  for 

then 

•KX.,  Y^)  = 

if  Y.  < X. 
otherwise 

j l<o  $‘ciion 

' C.:f  Ce.:;;on  □ 

r:: 

■ ! 
, - J 

1-  . 

u ■= 

m 

1 

i=l  j 

? 4>(X.  , Y.) 

-1  ^ ; 

1 

, L-‘  ' 

U'- 

is  the  well  known  two  sample  Mann-Whitney  statistic,  that  is 


U=  number  of  pairs  (X^ , Y^ ) such  that 


Y.  < X. 


Birnbaum  (1956)  showed  that  the  Mann-h'hi  tnry  statistic  U could  he 
used  to  estimate  1 - R (Probability  of  failure),  and  hence  R.  Tn 
particular 
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f = 1 - A *=  L/mn 


(2.1) 


was  proposed  as  an  estimator  of  P = Pr  (failure),  and  it  was  used 
to  obtain  one  sided  confidence  interval  for  P for  the  cases  F 
known,  G unknown  (m  ->■  ~) , and  both  F and  G unknown.  Birnbaum 
and  McCarty  (1958)  considered  a numerical  procedure  for  computing 
the  sample  sizes  needed  for  the  confidence  interval  based  on  U/mn. 

Owen,  Craswell  and  Hanson  (1964)  showed  that  the  assumption  of 
continuity  required  in  Birnbaum  (1956)  was  not  essential  and  produced 
some  tables  for  use  in  computing  sample  sized  and  confidence  intervals 
for  the  Birnbaum-McCarty  jirocedure . 

Govindarajulu  (1968)  also  has  explicitly  derived  one  sided  and 
two  sided  distribution  free  confidence  bounds  for  P based  on  the 
asymptotic  normality  of  ^ = U/mn.  This  bounds  are  approximately 
one  half  of  the  corresponding  bounds  due  to  Birnbaum  and  McCarty 
(1958).  In  particular,  Govinderajulu  showed  that  for  all  F and  G and 
large  m or  n,  the  solution  c of  the  equations 

P(P  s P + e)  = P(P  ^P-c)J:y»0<Y<1 
is  given  by 

-V2 

c ^ (4  V)  » (Y), 

and  the  solution  of  the  equation 

P(|P-P|se)^Y,  0<Y<1 

TOO 


V 


is  given  by 


c ^ 


(4  V) 


:^/2 


-1 

4> 


,1  + Y 


). 


Here 


4>(X)  = L , 

/2ir 

and  t is  the  inverse  function  of 

Recently  Govindara julu  (1974)  has  also  considered  a sequential  dis- 
tribution-free procedure  for  obtaining  fixed-width  confidence  limits 
for  P (and  hence  for  R) . However,  in  the  absence  of  additional 
numerical  computation,  it  is  not  known  how  good  is  the  performance 
of  this  sequential  procedure, 

3.  Bayesian  Approach 

Not  much  has  been  done  from  the  Bayesian  point  of  view 
Enis  and  Geisser  (1971)  investigated  Bayesian  approach  for  estimating 
R assuming  X and  Y to  be  independently  distributed  and  that  X and 
Y are  either  exponentially  distributed  or  normally  distributed. 

4 . Parametric  Approach 

In  many  situations,  the  distribution  of  X or  (of  both  X and  Y) 
will  be  known,  and  it  is  desired  to  obtain  parametric  solutions. 

Thus,  in  case  of  missile  flights,  the  stress  may  be  expensive  to 
sample,  but  the  physical  characteristics  of  the  missile  system,  such 
as  the  propulsive  force,  angle  of  elevation,  changes  in  atmospheric 
condition,  and  so  on  may  all  have  known  distributions;  consequently. 
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the  distribution  of  stresses  may  be  calculated.  In  this  section,  we 
shall  consider  the  problem  of  estimating  R (or  P)  for  specific  para- 
metric distributions. 


4.1  Normal  Distribution:  Owen,  Craswell  and  Hanson  (19C4)  considered 

above  problem  and  gave  one  sided  confidence  intervals  for  R when 
both  stress  and  strength  are  (a)  jointly  bivariate  normally  distribut- 
ed and  observations  are  in  pairs,  6r  (b)  when  X and  Y are  indepen- 
dent normal  with  a common  unknown  variance.  Note  if  X and  Y follow 
a joint  bivariate  distribution 


i 


where  y = EC^)  and  = E(Y  - Church  and  Harris  considered 

considered  the  estimator 


4>(V),  say, 


“ 1 ” 7 " - 7 

where  Y = - Z Y.  and  = v (y  . y)^/(n  - 1),  from  which  they 
i=l  ^ 1 1 

obtained  the  following  confidence  interval  for  R. 


P{$(V  - (j)‘^(l  - ^)  e^)  < R < $ (V  + $'^(1  - 5y)}=^l-T 


Similarly,  a one  sided  confidence  interval  is  given  by 


P{R  > t>  (y  - 4>‘^(1  - Y ) S^)}  = 1- Y . 


Here 


1 + 


( ^ - 

'■  n 


Y^s^ 


2(n  - IKl  + 5^)2 


The  confidence  interval  obtained  by  Church  and  Harris  compare 
favorably  with  that  of  Govindarajulu  (1968).  Their  procedure,  al- 
though empirically  demonstrated  to  be  superior  to  that  of  Govindarajulu 
is,  however,  inexact  since  it  uses  the  asymptotic  normal  approxima- 
tion of  a given  statistic  and  requires  the  substitution  of  the  popu- 
lation mean  and  standard  deviations  by  their  observed  sample  values. 

In  fact,  all  the  parametric  estimators  suffer  from  same  weakness  as 
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they  are  based  on  maximum  likelihood  estimators.  Mazumdar  (1970), 
has  considered  the  same  problem  of  obtaining  point  and  interval 
estimates  of  reliability  and  obtained  mvue  of  reliability  using 
interference  theory.  Minimum  variance  unbiased  estimator  of  R in 
the  normal  case  has  also  been  considered  by  Downton  (1973). 

4.2  Gamma  and  Exponential  distribution:  Since  in  many  physical 

situations,  specially  in  reliability  and  life  testing  problems,  ex- 
ponential and  gamma  distributions  provide  more  realistic  models, 
it  is  desirable  to  obtain  estimators  of  R in  these  cases. 

Let  X and  Y be  independently  distributed  with  density  func- 
tions 

f(x)  = 1——  xP'^  ' , X > 0,  p > 0 

r(p)aP 


g(y)  = ^ y^‘^  , y > 0,  q > 0 

F(q)  6^ 


respectively.  Then 


R P(X<Y)  = 


/ [l-G(x)]  dF(x) 

0 


/ 

0 


X r(q)B"^ 


e y*^  ^dy 


1 vP-1 


r(p)aP 


X*  dx 


k-0  r(p)r(k*l)  (a+r)P^^ 
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Here  p and  q are  assumed  to  be  Icnown  integers.  If  two  independent 

random  samples  and  , Y2 . • . from  the  two  gamma 

populations  are  available  mle  of  a and  3 are  given  by  a=^  and 
Y 

3=—.  Hence  mle  of  R is 

M 

j . ‘'j'  rcp*M  S'^gP  . 
k»o  r (p;r(kti) 

As  special  cases,  if  q=l,  that  is  if  X follows  the  gamma  distri- 
bution and  Y follows  the  exponential  distribution 

R = 

Finally,  if  both  p and  q are  equal  to  1 , we  have  the  case  of  two 
independent  exponential  distributions  and  we  have 


a + 8 X + Y 


The  distribution  of  for  large  m and  n,  can  be  shown  to  be 
normal  and  hence  asymptotic  confidence  interval  for  R can  be  ob- 
tained. 

•Tong  (1974,  1975)  has  obtained  mvube  of  R for  gamma  and  ex- 
ponential distributions.  The  variance  of  the  muvbe  of  R,  in  the 
exponential  case  has  been  derived  by  Kelley  et  al  (1976) 

4.3  V/cibui:..  distribution:  Let  X and  Y be  independent  random 

variables  each  followii g the  Weibull  distribution  with  common  shape 
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, a > 0 , X > 0 
, 3 > 0 , y > 0, 

We  can  readily  see 

R ■=  PCX"^  < = P(X  < Y)  = — 

a + 3 

Note  above  is  independent  of  6.  Again,- we  can  obtain  the  mle 
of  R to  be 

R = 3/(Q  + 3) 

whore  a and  3 are  mle  of  a and  3- 

4.4  Bivariate  exponential  distribution:  Since  exponential  distri- 

bution is  considered  a useful  model  in  life  testing  problems,  it  is 
desirable  to  consider  bivariate  analogue  of  univariate  exponential 
distributions  which  will  have  properties  similar  to  the  univariate 
exponential  distribution.  Marshall  and  Olkin  (1967)  have  proposed 
a very  important  bivariate  exponential  distribution  (BVE) , which 
is  given  by 

F(x,y)=PCX>^c,Y>y)-e'^l’''^2^'^lf"^*^^*>'\o<X^,X2,X^2^»,Xj+Aj2>0,A2+A^2=^0(x>0,y>0). 


parameter  6,  That  is  let 


F(x)  = 1 - e 


-X 


fa 


6/ 

G(y)  = 1 - e‘^  /3 
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The  BVh  does  arise  in  several  natural  ways  and  its  properties  ap- 
pear to  be  fundamental.  In  particular,  marginal  distributions  of 
BVI;  are  exponential  and  BVE  has  the  loss  of  memory  property  (BMP) 
given  by 


F(x  + t,y  + t)  = F(x,y)F(t,t)  for  Sj^,  S2>  tsO 


However,  this  distribution  is  not  absolutely  continuous  and  there 
are  clearly  situations  when  it  can  not  be  applied.  Thus,  from  data, 
it  is  found  that  X^Y  for  any  pair  (X,Y)  the  model  is  clearly 
not  applicable.  An  alternative  absolutely  continuous  distribution 
related  to  the  BVE  and  having  some  of  its  properties  would  appear 
to  be  of  interest.  To  this  end.  Block  and  Basu  (1974)  have  proposed 
an  absolutely  continuous  bivariate  exponential  extension  (ACBVE) , 
which  turns  out  to  be  the  absolutely  continuous  part  of  the  BVE  of 
Marshall  and  Olkin.  ACBVE  is  also  seen  to  be  a variant  of  the 
distribution  Freund  (1961) . The  ACBVE  is  given  by 


F(x,y)  = 


X1  + X2 


X1  + X2 


exp  [ - Xj^x- X2y  - X (x ,y)  ] 


exp [ - Xmax (x ,y) ] for  x>0,y>0, 


Here 


X - Xf  + X2  + Xj2 


Estimates  of  R when  the  underlying  distribution  is  EVE  or 
ACBVE  has  been  obtained  by  Basu  (1976) . These  results  will  be  com- 


municated elsewhere. 

5 . Reliability  of  complex  systems 

The  model  described  before  can  be  extended  to  more  complex  sys- 
tems. For  example,  a single  component  system  of  strength  Y could 
be  subjected  to  k different  independent  stresses  ,X2  , . . . Xj^ . 

Here  reliability  of  the  system  is  given  by 

R = P{Xj<Y,  X2<Y,  . . . ,Xj^<Y} 
or 

R = P{max(X^,X2 Xj^)  < Y}. 

f 

I 

l 

An  example  of  interest  is  the  case  where  a beam  with  strength  Y is 

subjected  to  several  stresses  Xj^,X2 Xj^.  Another  similar  pro- 

: blem  of  interest  is  to  evaluate  the  reliability  function  R'  of  a 

k -component  system  of  strengths  •••  » respectively 

each  of  which  is  subject  to  a common  stress  X.  Here 

R = P{X<Yj,  X<Y2,  ...  , X<Yj^} 

= P{X<min  (Yj^ ....  , Yj^) } . 

As  an  example,  the  flow  of  a current  X through  an  electronic  com- 
ponent assembled  from  several  subcomponents  with  abilities  to  accom- 
modate currents  X^,Y2 Yj^  would  follow  this  pattern. 
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Chandra  (1975)  has  considered  the  problem  of  estimating  R and 
R'  under  the  assumption  that  the  X's  and  Y's  are  all  independent 


random  variables  and  (a)  all  follow  normal  distributions,  (b)  Y's 
are  all  exponential  and  X is  normal  with  known  variance. 

Bhattacharyya  and  Johnson  (1974)  considered  the  problem  of 
estimating  reliability  function  R for  a more  complex  m-out-of-k 
system.  Here  each  of  m components  of  a system  of  strengths  Yj^ , 

Y2 Yj^  is  subjected  to  a stress  X and  the  system  survives  if 

at  least  m out  of  the  k components  survive.  Assuming  X,Yj^,  ...  ,Yj^ 
to  be  independent  with  distribution  functions  F(x),  Gj^fy^^),  02(72)1 
...  iGj^(7jj).  Bhattacharyya  and  Johnson  considered  the  problem  of 
estimating  the  reliability  function  R=P^(at  least  m of  the  Y^^ , .... 
Yj^  exceed  X),  under  the  assumption  Gj^  = G2= . . =Gj^=G,  say>  and  that  F 
and  G are  exponential  distributions  with  known  scale  parameters. 

Here 

R=  E M /“  [l-G(x)HG(x)]^'“d(F(x))  . 
a=m  ' ' 

Bhattacharyya  and  Johnson  (1973)  have  also  considered  a nonparametric 
approach  for  the  above  problem. 

The  author  is  currently  investigating  additional  problems  in 
this  area  results  of  which  will  be  communicated  elsewhere. 
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